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"-•*  Open  cavities  on  aircraft  exposed  to  high  speed  flow,  such  as  bomb 
bays  can  give  rise  to  intense  self-induced  pressure  oscillations.  The 
amplitude  of  these  oscillations,  under  certain  flight  conditions,  can 
cause  structural  damage.  Substantial  experimental  and  analytical  efforts 
have  investigated  these  pressure  fluctuations,  resulting  in  some  under- 
standing of  the  complex  interaction  of  the  external  shear  layer  and 
cavity  acoustical  disturbances.  However , no  numerical  computations  have 
been  obtained  for  tine  complete  governing  fluid  mechanical  equations.  The 
purpose  of  this  study  is  to  obtain  numerical  solutions  of  the  Navier-Stokes 
equations  for  an  open  cavity  in  order  to  provide  a now  tool  for  the  analysis 
of  this  phenomenon. 
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NOMENCLATURE 


speed  of  sound 

coefficient  in  pressure  perturbation  equation 

complex  propagation  speed 

specific  heat  at  constant  pressure 

cavity  depth 

specific  internal  energy 

vector  fluxes 

frequency  of  wave 

c^/U^  propagation  velocity  ratio 

cavity  length 

mode  number 

Mach  number 

node  number 

pressure 

heat  transfer  rate 
gas  constant 
Reynolds  number 
time 

temperature 

velocity  components  in  Cartesian  frame 

vector  of  dependent  variables 

Cartesian  coordinates 

2nd /A  dimensionless  wave  number 

ratio  of  specific  heats 

shear  layer  thickness 
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X 

V 

p 

0 

‘i 

4 

01 


wave  length 

viscosity  (molecular  and  eddy) 

density 

normal  stress 

viscous  shear  stress 

amplitude  of  perturvation  velocity 

2n£  frequency 


Subscripts 


00 

0 
W 

r 

1 
1 

2 


frcestream  condition 
stagnation  condition 
wall  condition 
rei>l  part 
imaginary  part 
forward  traveling  wave 
rearward  traveling  wave 


Superscripts 

• instantancuos  perturbation  variable 

- - vector 
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I.  Background 


As  early  as  1955,  Krishnamurty  investigated  flow  induced  pressure 

oscillations  in  open  cavities.  Other  investigators^ »3>^>5»6>7,8>9  have 

conducted  extensive  research  in  an  attempt  to  understand  the  physical 

3 

mechanisms.  Heller  and  Bliss  used  a water  table  to  simulate  supersonic 
airflow  over  open  cavities.  They  found  that  the  inherently  unstable 
shear  layer  fluctuates,  causing  periodic  mass  addition  and  expulsion  from 
the  cavity  (Figure  1) . When  the  rear  reattachment  point  of  the  shear 
layer  enters  the  cavity,  a stagnation  point  is  created;  thus  increasing 
the  local  cavity  pressure.  This  mass  addition  creates  a traveling  pres- 
sure wave  (as  in  a shock  tube) , which  moves  forward  in  the  cavity  (at 
supersonic  speed  relative  to  free  stream),  trailing  an  oblique  shock  in 
the  free  stream.  When  the  traveling  shock  wave  reflects  from  the  forward 
bulkhead,  a pressure  doubling  occurs  in  the  cavity  while  disturbances  in 
the  external  flow  are  not  reflected,  and  thus  a pressure  jump  across  the 
shear  layer  deflects  the  shear  layer.  The  reflected  traveling  shock  wave 
in  the  cavity  is  now  moving  at  subsonic  speed  relative  to  the  free  stream, 
hence,  generates  no  oblique  shock  wave  in  the  free  stream.  As  die  cavity 
traveling  shock  wave  approaches  the  rear  bulkhead,  the  shear  layer  bulges 
outward,  and  mass  is  ejoctod  out  of  the  cavity.  The  entire  process  then 
repeats  itself  in  a periodic  fashion. 

Thus,  sufficient  experience  from  such  extensive  measurements  exists 
so  that  a qualitative  description  of  the  flow  process  can  be  obtained. 
However,  a quantitative  prediction  method  does  not  exist  which  is  the 
motivation  for  the  present  investigation. 
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To  study  self-induced  pressure  oscillations  in  an  open  cavity,  an 
analytic  study  was  first  accomplished,  followed  by  a numerical  computa- 
tion of  the  Navier-Stokes  equations  and  a comparison  with  previous  experi- 
mental investigations. 

II.  Analytical  Study 

Consider  a system  of  traveling  waves  which  produce  the  wave-diagram 
(x,t)  shown  in  Figure  2.  A resonant  situation  arises  when  a forcing 
function  excites  the  shear  layer  in  the  frequency  range  where  amplifica- 
tion is  possible.  The  disturbances  will  grow  until  a limit  xycle  is 
reached  due  to  viscous  dissipation.  A standing  wave  exists  in  the  cavity 
when  both  the  upstream  and  downstream  traveling  waves  are  synchronized. 
This  wave  pattern  of  Figure  2 may  be  approximated  by  considering  forward 
and  rearward  traveling  pressure  waves  of  equal  intensity  but  different 
propagation  velocities  and  wave  numbers. 


ia.(x  - c.t)  ia,(-x  - c.t) 

p “ Ae  + Ae 


U) 


The  frequency  of  the  pressure  pulse  can  be  determined  directly  from  the 
wave  diagram. 


a.i  a 

f c4  c2 


(2) 


where  m is  the  number  of  waves  or  mode  number. 

From  observation  of  cavity  oscillations  the  rearward  traveling  wave 
(cig)  is  known  to  be  att  acoustical  disturbance  traveling  at  the  speed  of 
sound  in  the  cavity. 

c,  - Vm  A ♦ 0.2N* 
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The  progagation  speed  of  the  forward  traveling  wave  has  been  observed  to 
be  about  half  of  the  free  stream  valued 


1/2 


The  frequency  can  be  evaulated  by  using  these  results. 

mU 

00 

f a Rossiter  s Formula  (Ref  4) 

L(Mq  + k'1) 


(A) 


(5) 


This  equation  have  been  used  successfully  for  determining  cavity  resonant 
frequencies  but  no  satisfactory  prediction  method  has  been  available  for 
determining  the  disturbance  intensity  of  the  different  inodes. 

The  mode  shape  of  these  standing  waves  can  also  be  deduced  from  the 
preceding  equations.  For  a standing  wave  to  occur  both  waves  must  possess 
the  same  frequency. 


tt  **  al  C1  * a2  c2  " 2,Tf 


(6) 


Combining  this  result  with  equation  2 produces  the  following  relationship: 

al  + a2  * ~ (7) 


The  mode  shape  may  be  obtained  by  utilizing  this  information  in  the  pres- 
sure equation  (equation  1)  and  computing  the  rtas  value  of  pressure  over  a 
complete  cycle. 


p 

rrms 

p 

rrms 


A(i  + cos  (a^  + a,)*) * 
mnxl 


cos 


(8) 


These  patterns  have  been  documented  in  Ref  3.  Nodes  will  occur  when 


m*x  tt 

— * n2{ 


n " odd 
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Mode  m 

1 

Nodes  n 
3 5 

1 

1/2 



2 

1/A 

3/A 

3 

1/6 

3/6  5/6 

Table  I Node  Location  for  Various  Modes 
HI.  Stability  Analysis 

A necessary  condition  for  resonance  is  that  one  of  the  waves  must  be 

unstable  for  the  oscillation  to  persist,  otherwise  the  disturbance  will 

dissipate  after  an  initial  transient.  Mathematically  this  means  that  the 

wave  speed  is  complex,  i.e.,  c ■ c^  + ici  with  c^>  0 unstable.  The 

stability  of  the  shear  layer  will  now  be  examined. 

Rayleigh*^,  in  1880,  showed  for  inviscid  incompressible  flow  that 

velocity  profiles  with  inflection  points  are  unstable.  Recently,  Michalke** 

confirmed  that  a shear  layer  is  unstable  but  only  at  low  frequencies 

(A/6  > 2n  or  fGAtt/tJ^  < 1).  It  was  felt  that  more  information  about  the 

stability  of  a compressible  shear  layer  was  needed,  therefore,  a linear 

12 

stability  analysis  was  undertaken  . the  governing  Euler  equations  were 
linearized  by  assuming  small  perturbations  caused  by  small  amplitude 
traveling  waves.  The  resulting  stability  equation  first  derived  by  Lees 
and  Lin^,  reduces  to  the  Rayleigh  equation  for  incompressible  flow. 

•2(g~l(U  - c)$y  - g"1  Uy  $)y  - a2(U  - c>$  (9) 
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where 


v - <Ky)eia<x  " ct)  (10) 

g “ a2  - (U  - c)2  (11) 

12 

The  eigenvalues  of  this  Rayleigh  equation  were  then  computed  for 
a shear  layer  with  a hyperbolic  tangent  velocity  profile. 


jj—  «*  0.5  (1  + tanh  |-) 


(12) 


The  propagation  velocity  of  the  disturbances  is  shown  in  Figure  3. 


c 

r 


= k 


(13) 


This  is  the  k value  determined  experimentally  by  Rossiter^  and  found  to 
be  in  excellent  agreement  with  his  results. 


.5  < k < .6 

The  amplification  factors  (c^  > 0)  were  found  to  be  a function  of  wave 
number  (a)  and  Mach  number  (Figure  4).  Instability  was  observed  only  i\  r 
wave  numbers  less  than  unity.  This  implies  short  cavities  (L  < 2nd)  will 
not  resonate.  Note  the  Rayleigh  instability  vanishes  above  M * 2.5.  This 
result  confirms  previous  experimental  and  numerical  results  that  separated 
flown  are  more  stable  at  supersonic  speeds  titan  at  subsonic. 

It  is  possible  to  predict  the  relative  intensity  of  the  different 
modes  occurring  in  an  open  cavity.  Consider  a cavity  of  L * 91.44  cm  and 
* 0.85,  l)w  * 286.$  mps 

Hence, 

—L 

y * k * .52  (from  Fig  3) 

m 

M - 0.79 

f • 115  (lit)  Rossi  tor's  Eqo. 
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These  different  mode  frequencies  are  plotted  in  Figure  5 using 
results  of  Ref  12,  and  the  second  mode  is  observed  to  have  the  greatest 
amplification  while  modes  four  and  above  are  found  to  be  stable.  A 
spectral  analysis  of  a wind  tunnel  test  (Ref  1A)  for  an  open  cavity  at 
these  same  conditions  is  also  shown  in  Figure  5.  Observe  that  only  the 
first  four  modes  are  dominant  and  the  relative  amplitudes  of  these  four 
inodes  are  consistent  with  the  analytic  amplification  factors. 

The  intensity  of  the  pressure  fluctuation  will  be  proportional  to 
the  following: 

ac 

P™c  W - q(e  - 1)  (1A) 

nas 

aci 

Selecting  the  peak  value  of  e for  each  Mach  number  and  multiplying  by 

q/P  (M)  the  relative  intensity  as  a function  of  Mach  number  may  be  deduced 
o 

(Fig  6).  The  peak  pressure  value  occurring  in  a series  of  wind  tunnel 
tests  at  different  Mach  numbers  can  be  expected  to  occur  near  Mach  one. 
This  is  confirmed  in  Reference  3 and  1A. 

XV.  Summary  of  Analytic  Results 

The  analytic  results  based  primarily  upon  stability  theory  provide 
us  with  the  following  conclusions. 

a.  Shear  layers  (with  inflection  points  in  the  velocity  profile) 

•re  unstable  but  only  for  low  frequencies;  fd/l^  < 1/An. 

b.  Short  cavities  (L  < 2«S)  will  not  resonate. 

c.  No  Rayleigh  instability  occurs  above  Mach  number  2.5. 

d.  Peak  amplification  Occurs  at  about  half  the  cut-off  frequency 
creating  a situation  where  modes  other  than  the  fundamental  can  dominate. 


2 16 


e,  Maximum  pressure  intensity  of  a shear  layer  oscillation  will 
occur  near  Mach  1 (in  a wind  tunnel  with  constant  pQ) . 

Linear  stability  theory  is  therefore  useful  in  explaining  the  cause 
of  the  resonance  and  in  estimating  the  relative  intensity  of  the  various 
modes  including  the  influence  of  Mach  number.  However,  the  full  non- 
linear equations  are  required  to  determine  the  absolute  level  of  the  pres- 
sure intensity.  For  that  reason  the  numerical  solution  of  the  unsteady 
Navier-Stokes  equations  will  be  considered  next. 

V.  Numerical  Computation 

With  the  completion  of  the  simplified  analytic  approach  (ie,  inviscid, 
linear  stability  theory)  a numerical  solution  of  the  exact  equations  was 
attempted  to  further  improve  Che  prediction  capability.  The  analysis 
served  to  identify  the  primary  mechanism  involved  in  the  oscillation  and 
greatly  assisted  in  the  determination  of  the  grid  point  distribution  and 
step  sizes  required  to  resolve  the  flow  features. 

A case  to  compute  was  selected  for  which  experimental  data  were 
available.  Thu  test  conditions  of  Heller  and  Bliss'*  were  selected  to 
compare  the  numerical  computations  (Figure  7). 

H - 1.5  L ■ 91.44  cm 

L/D  - 2.25  D * 40.64  cm 

Re  « 1.28  x 10^/b  6 * 2.54  cm 

o 

The  width  of  the  cavity  was  22.66  cm  and  found  not  to  be  a major 
factor  in  the  overall  phenomenon.  For  this  reason,  a two-dimensional 
computation  appeared  to  be  justified  for  the  initial  studies. 
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VI.  Governing  Equations 


The  time  dependent  explicit  finite  difference  method  originated  by 
R.  MacCormack1^  was  selected  to  perform  the  numerical  calculations.  The 
two-dimensional  Navier-Stokes  equations  follow: 


3U 

3t 


3E 

+ 35  + 
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3y 


(15) 
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where 


°xx  " “P  ~ WvN  * u + 2y  u^ 


t * u(u  v ) 
xy  y x 


oyy  - -p  - 2/3pV  • u 2puy 


(16) 


The  turbulent  closure  of  the  present  problem  was  achieved  by  implement- 
ing the  Cebcci-Smith  eddy  viscosity  model  with  relaxation  modification16. 

The  relaxation  turbulence  model  was  used  in  an  attempt  to  describe  the 
adjustment  of  the  turbulence  structure  from  an  attached  boundary  layer  to 
an  oscillatory  free  shear  layer.  The  relaxation  length  scale  was  assigned 
a value  of  55  boundary  layer  thicknesses.  Since  there  is  no  guidance  to 
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assess  the  accuracy  of  the  turbulence  model  relative  to  low  frequency 

fluctuations,  a parametric  study  seemed  to  be  necessary.  In  the  present 

analysis,  several  consecutive  calculations  with  suppressed  eddy  viscosity 

were  performed  and  the  numerical  results  exhibited  only  a minor  departure 

from  the  basic  solution.  Hence,  the  present  eddy  viscosity  model  was  felt 

16  17 

to  be  adequate.  The  computer  program  previously  used  by  Shang  ’ was 
modified  to  include  the  appropriate  boundary  conditions  for  this  problem. 


VII.  Boundary  Conditions 

Four  faces  require  attention  in  the  specification  of  boundary  con- 
ditions (Figure  7). 


Wall  and  Cavity  Surfaces: 


On  solid  surfaces,  the  velocity  components  vanish,  and  the  wall 
temperature  must  be  prescribed.  In  addition,  the  pressure  is  derived 
from  the  respective  compatibility  conditions  of  the  Koraent-.na  equations. 


u * 0 


v « C 


T « T 
w o 


I*  " lx  “ |y  (w(uy  + V.)) 


(IS) 


§y.-“  l-2/3'^  * u + *»»y)  - n lw(uy  + V1 


Upstream  Condition: 

A supersonic  free  stream  with  a known  boundary  layer  profile  is  given. 

2 

u * u (y)  T - T - 


2C 


P 


(19) 


v * 0 


P ■ P, 


w 
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Downstream  Condition: 


A mild  boundary  condition  is  prescribed  to  avoid  major  reflections 
of  disturbances. 


0 


A similar  no-reflection  condition  is  adopted. 


Where  is  the  outgoing  characteristics  on  the  upper  boundary  of  the 
computational  domain. 


Initial  Condition: 

The  upstream  condition  is  imposed  as  the  initial  condition  for  the 
flow  outside  the  cavity.  Inside  the  cavity  initially  the  flow  is  assumed 
to  be  static. 


u « 0 


v » 0 T ■ T 


P ■ P« 


VIII.  Numerical  Procedure 


i 1 

MacCormack'a  J alternating-direction-explicit  numerical  scheme  was 
adopted  for  the  present  analysis.  For  this  case  pressure  damping  was 
required  due  to  the  transient  multi-wave  structure  occurring  in  the  flow 
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field.  The  current  philosophy  in  computational  ftuid  dynamics  is  to 
employ  a body  oriented  coordinate  system  which  turns  out  to  be  Cartesian 
in  this  case  with  non-uniform  step  si2e.  For  cases  presented  here,  a 
grid  of  78  x 52  was  used  to  represent  a field  size  of  182.9  cm  x 91.44  cm 
enclosing  a 91.44  cm  x 40.64  cm  cavity  (see  Figure  7).  In  order  to 
achieve  the  desired  temporal  resolution,  a time  step  corresponding  to  a 
Courant  number  of  0.2  was  used.  All  calculations  were  performed  on  a CDC 
6600  computer.  The  data  processing  rate  was  0.0017  sec  per  grid  point 
per  time  step.  The  central  core  memory  required  for  the  present  problem 
is  205K  octal. 

Although  dispersion,  dissipation  and  phase  errors  are  not  negligible 

with  the  step  sizes  employed,  previous  numerical  investigations  of  viscous 

interaction  problems”  ^ using  comparable  step  sizes  have  shown  good 

agreement  (+5  widr  experimental  data  for  the  most  significant  features 

of  the  flow.  In  particular,  periodic  motions  around  a transonic  airfoil 

18 

have  been  studied  by  Levy  with  a basic  MacCormack's  scheme.  His  results 
exhibited  good  agreement  with  experimental  data  not  only  in  the  pattern  of 
a simple  wave  train  but  also  in  the  predicted  reduced  frequency.  Therefore, 
no  additional  modification  other  than  a simple  controlled  spatial  averaging 
was  used  to  correct  the  possible  dispersion  error  for  the  compound  wave 
problem  investigated. 


IX.  Discussion  of  Results 

Time  dependent  numerical  computations  of  supersonic  flow  over  an 
open  cavity  were  accomplished  utilizing  MacCormack's  finite  difference 
explicit  method.  The  entire  velocity  field  over  the  cavity  is  shown  in 
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Figure  8.  Since  the  problem  is  strictly  a time  dependent  phenomenon,  only 
a typical  velocity  distribution  is  presented  here  (t  = 0.0062  sec)  to 
reveal  the  basic  features.  The  most  obvious  feature  is  that  the  flow 
field  within  the  confined  cavity  is  subsonic,  except  perhaps  the  region 
adjacent  to  the  cavity  opening.  The  experimental  investigation  also 
recorded  the  identical  observation.  The  orderly  development  of  the  shear 
layer  above  the  cavity  is  also  clearly  exhibited.  An  attached  turbulent 
boundary  layer  upstream  of  the  cavity  separates  at  the  forward  bulkhead 
to  form  a free  shear  layer  over  the  cavity  and  finally  reattaches  down- 
stream of  the  cavity.  Due  to  the  smaller  magnitude  of  the  velocity  com- 
ponents within  the  cavity,  the  velocity  distribution  could  not  be  shown 
with  the  same  scale  as  that  of  the  outer  shear  layer.  A magnified 
velocity  profile  in  the  cavity  is  presented  in  Figure  9.  All  velocity 
distributions  were  drawn  at  a scale  ccn  times  greater  than  that  in  Figure 
8.  A recirculation  flow  configuration  is  demonstrated  with  the  center  of 
the  recirculation  located  near  the  upper  corner  of  the  rear  bulkhead.  For 
clarity  only  every  fourth  velocity  point  in  the  stveamwise  direction  was 
presented  in  both  figures  (Fig  8 and  9). 

A quantitative  comparison  of  the  calculated  mean  velocity  distribution 
with  experiment  is  presented  in  Figure  10.  The  calculated  velocity  pro- 
file immediately  upstream  and  downstream  of  the  cavity  compares  very  well 
with  the  experimental  measurements.  The  upstream  velocity  profile  is 
essentially  independent  of  time  as  expected,  and  nearly  duplicates  the 
data.  The  downstream  velocity  profile  indicates  that  the  reattached  shear 
layer  thickens  significantly  over  the  cavity.  The  difference  between  data 
and  calculation  is  a mere  seven  percent.  The  computed  velocity  profile 
downstream  of  the  cavity  exhibits  an  oscillatory  behavior  in  the  inner  por- 
tion of  the  boundary  layer  which  is  confirmed  by  the  experimental  observation 
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A comparison  of  the  time-average  surface  pressure  distribution  in  the 

cavity  was  performed.  In  Figure  11,  one  observes  that  the  computed  mean 

pressure,  normalized  by  the  free  stagnation  pressure,  uniformly  under  pre- 

3 

diets  the  experimental  data  by  about  10  percent.  The  maximum  deviation 
between  data  and  calculation  occurs  near  the  rear  bulkhead  were  the  pres- 
sure difference  accross  the  cavity  opening  and  cavity  floor  also  reaches 

a maximum.  Nevertheless,  the  calculated  results  indicate  the  identical 

3 

trend  to  that  of  the  experiment  . 

In  Figure  12,  a history  of  the  static  pressure  at  x/L  = 0.66  and 
y/L  **  0.960  is  monitored  (y/L  = 0.960,  x/L  = 0.33,  0.50  and  0.66).  Per- 
sistent oscillatory  static  pressures  appeared  within  one  characteristic 
time,  t The  characteristic  time  is  defined  as  the  period  of  time 
required  for  a fluid  particle  to  traverse  the  length  of  the  cavity  at 

_3 

freestreara  speed.  Tor  the  present  problem  t^  has  a value  of  1.987  x 10 
seconds.  Since  the  pressure  oscillation  over  the  cavity  is  composed  of 
several  frequencies  of  diff event  amplitude,  a compound  wave  system 
develops.  The  present  calculation  was  carried  out  only  over  a time  span 
of  about  nine  characteristic  times  (t  «=  1.82  x 10  secs).  This  result  is 
compared  with  oscillogram  data  for  the  test  results  of  Ref.  IS  in  Figure 
12.  In  order  to  permit  a qualitative  comparison,  computed  results  aie 
repeated  for  several  cycles.  One  can  decert  certain  similarity  between 
the  data  and  the  present  result.  In  prii.*  pie,  the  solid  surface  con- 
straint and  wave  interference  phenomenon  r*  che  present  problem  are 
Inherently  nonlinear.  Therefore,  raut_on  must  be  exercised  in  discerning 
the  discrete  frequencies  between  the  fundamental  tu.des  of  oscillation. 
Other  basic  information  such  as  the  rclutivc  phase  angle  aud  the  amplitude 
of  each  distinctive  wave  requires  further  spectral  analysis. 
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A spectral  analysis  of  this  compound  wave  probably  is  the  only 
reliable  means  for  accurately  obtaining  mode  frequency.  However,  this 
is  impractical  due  to  the  large  amount  of  computer  time  required  to 
obtain  solutions  for  a sufficiently  long  duration.  Hence,  an  analysis 
was  accomplished  by  assuming  the  waves  to  be  commensurable.  The  spectral 
analysis  reveals  four  distinctive  discrete  frequencies  of  154  Hz,  308  Hz, 
462  Hz,  and  616  Hz  recognized  as  the  first,  second,  third  and  fourth  mode 
respectively  of  the  oscillatory  pressure  disturbance.  The  higher  modes 
of  oscillation  decay  rapidly  as  one  may  observe  in  Figure  13.  Good 
agreement  between  the  experimental  measurement  and  present  result  is 
observed.  Both  exhibit  a dominate  second  mode  of  the  pressure  oscillation. 
The  level  of  pressure  oscillation  in  db  can  be  evaluated  as 

p 

p(db)  » 20  log  + 189 

%> 

where  is  the  dynamic  pressure  (54  K pascal).  The  detected  frequency 

for  the  second  mode  (308)  compares  well  with  the  experimental  data  (300  Hz) 

and  Rossiter's  prediction^  (328  Hz).  The  fluctuating  pressure  level 

between  the  data  and  present  result  is  within  about  ton  percent. 

The  compound  wave  pattern  is  best  illustrated  in  Figure  14.  The 

propagation  of  the  wave  train  from  the  forward  bulkhead  is  presented  for 

-3 

a fixed  time  interval  of  0.64  x 10  seconds.  One  observes  the  rearward 
traveling  propagation  wave  has  an  unmodulated  amplitude  until  interacting 
f with  reflected  waves  from  the  rear  bulkhead.  No  repeatable  wave  front  can 

j be  identified  downstream  of  x/L  **  0.75.  Two  pieces  of  important  informa- 

i 

tion  have  been  determined  from  this  graph,  namely  the  rearward  wave  pro- 
pagation speed  and  the  amplitude  of  the  pressure  oscillation.  The  pre- 
dicted wave  speed  has  a value  of  244.4  mps  (or  k ■ .53)  and  is  in  agreement 

I 

V: 
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with  the  predicted  value  from  Figure  3.  The  amplitude  of  the  oscillating 
pressure  also  agrees  well  with  data,  (+  6.464  K pascal  vs  7.182  K pascal; 
or  170  db)  with  the  discrepancy  about  10%. 

In  Figure  15  the  sequence  of  density  contours  from  the  numerical  com- 
putation is  shown  for  a complete  cycle  of  the  periodic  motion.  The  for- 
ward and  rearward  moving  wave  system  originated  from  the  instability  of  the 
free  shear  layer  and  the  reflection  at  the  rear  bulkhead  can  be  easily 
recognized  in  the  cavity.  The  generation  and  movement  of  the  external 
shock  wave  system  also  can  be  recognized.  These  compare  favorably  with 
the  wave  pattern  shown  in  Figure  1 for  the  water  table  experiment. 

X.  Conclusions 

The  pressure  oscillation  for  supersonic  flow  over  an  open  cavity  has 
been  predicted  by  numerically  solving  the  unsteady  Navier-Stokes  equations. 
Both  the  predicted  frequency  and  magnitude  of  the  unsteady  pressure 
fluctuations  were  qualitatively  confirmed  through  experiment.  However,  a 
spectral  analysis  of  a numerical  solution  of  longer  duration  is  required 
for  complete  verification.  This  is  the  first  time  a complete  viscous 
solution  of  the  pressure  oscillating  cavity  tuis  been  obtained  and  displays 
the  outstanding  capability  inherent  in  the  numerical  methods  of  today. 
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Figure  4.  Amplification  Factor  v.s.  Wave  Number 
for  Different  Mach  Number  (Ref,  12) 
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Figure  5.  Comparison  of  Amplification  for  Different 
Frequencies  with  Experimental  Data 
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Figure  6.  Relative  Intensity  of  Disturbances  for  Different  Nach  Numbers 
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Figure  10.  Comparison  o £ Velocity  Profiles 
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Figure  13.  Comparison  of  Spectral  Analysis  Results 


